Graviton production in anti-de Sitter braneworld cosmology: 
A fully consistent treatment of the boundary condition 
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In recent work by two of us, [Durrer & Ruser, PRL 99, 071601 (2007); Ruser & Durrer PRD 
76, 104014 (2007)], graviton production due to a moving spacetime boundary (braneworld) in a 
five dimensional bulk has been considered. In the same way as the presence of a conducting plate 
modifies the electromagnetic vacuum, the presence of a brane modifies the graviton vacuum. As 
the brane moves, the time dependence of the resulting boundary condition leads to particle creation 
via the so called 'dynamical Casimir effect'. In our previous work a term in the boundary condition 
which is linear in the brane velocity has been neglected. In this work we develop a new approach 
which overcomes this approximation. We show that the previous results are not modified if the 
brane velocity is low. 
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I. INTRODUCTION 

The idea that our Universe is a 3 + 1 dimensional 
membrane in a higher dimensional 'bulk' spacetime has 
opened new exciting prospects for cosmology, for reviews 
see [H, 01- In the simplest braneworlds motivated by 
string theory, the standard model particles are confined 
to the brane and only the graviton can propagate in the 
bulk. Of particular interest is the Randall- Sundrum (RS) 
model , where the bulk is 5-dimensional anti-dc Sit- 
ter space. If the so called RS fine tuning condition is sat- 
isfied, it can be shown that gravity on the brane 'looks 
4-dimensional' at low energies. 

Within this model, cosmological evolution can be inter- 
preted as the motion of the physical brane, i.e. our Uni- 
verse, through the 5d bulk, acting as a moving boundary 
for bulk fields, in particular for 5d gravitational perturba- 
tions. Such a time-dependent boundary does in general 
lead to particle production via the dynamical Casimir 
effect [1,11]. 

Of course one can always choose coordinates with re- 
spect to which the brane is at rest, e.g. Gaussian nor- 
mal coordinates. This leads to a time dependent bulk 
resulting in the same effect, particle production from 
vacuum due to a time varying background metric. But 
then, usually (except in the case of de Sitter expansion 
on the brane [7[), the perturbation equation describing 
the evolution of gravitons is not separable and can only 
be treated with numerical simulations [1, 0, [l(| . Further- 
more, in a time dependent bulk, a mode decomposition 
is in general ambiguous and one cannot split the field in 
a zero mode and Kaluza-Klein (KK) modes in a unique 
way. One of the advantages of the dynamical Casimir 
effect approach presented in [lj], [l2j] is that it allows for 
a clear physical interpretation and in addition exhibits 
an analogy with quantum electrodynamics. 

Based on the picture of a moving brane in AdSs, we 
have studied graviton production in an ckpyrotic type 



scenario [13| where our Universe first approaches a sec- 
ond static brane. After a 'collision' the physical brane 
reverses direction and moves away from the static brane, 
see Fig. [T] For an observer on the brane, the first phase 
corresponds to a contracting Universe, the collision rep- 
resents the 'Big Bang' after which the Universe starts 
expanding (see Fig. 1). We do not model the details of 
this collision, but assume that the brane distance is still 
finite at the collision. This corresponds to a cutoff of 
all the physics which happens at scales smaller than the 
minimal brane distance when contraction reverses into 
expansion. In our results we assume this to be of the 
order of the string scale. We cut off the spectra at the 
string scale. This is a conservative assumption which sig- 
nifies that we neglect all the particle creation at energies 
higher than this scale. 

We have obtained the following important results in 
our previous papers [111, [12] : first of all, the energy den- 
sity of KK gravitons in AdSs scales like stiff matter, 
oc a~ 6 , where a denotes the scale factor introduced in 
Eq. ([2]). Therefore, KK gravitons in AdSs cannot rep- 
resent the dark matter in the Universe 1 . We have also 
seen that in the early Universe the back reaction from KK 
gravitons on the bulk geometry is likely to be important. 
Finally, we have derived a limit for the maximal brane 
velocity, the bounce velocity, w max <, 0.2 in order not to 
over-produce zero-mode (i.e. 4d) gravitons, the energy 
density of which is constrained by the nucleosynthesis 
bound. We have also calculated the spectra of both, the 
zero-mode and the KK gravitons. 

In this previous work we have, however, neglected 
a term linear in the brane velocity v in the boundary 
conditions (junction conditions) for the tensor pertur- 
bations. Here we derive a method which includes this 



1 See 1 1 4| for a discussion on a contradicting result in the literature. 
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FIG. 1: Two branes in an AdSs spacetime. The physical 
brane, a Friedmann universe with energy density p is on the 
left. While it is approaching the static brane its scale factor 
is decreasing, the Universe is contracting, and when it moves 
away from the static brane the Universe is expanding. L is the 
AdS curvature radius which is related to the brane tension T 
via Eq. |SJ|. The value of the scale factor of the brane metric 
as a function of the extra dimension y is also indicated. 



term and allows to treat the problem without any low 
velocity approximation. We show that the low velocity 
results previously obtained are not modified. Especially, 
the nucleosynthesis bound on the maximal brane veloc- 
ity, fmax~0.2, remains valid. In a subsequent study we 
shall investigate graviton production from branes which 
achieve high velocities in detail [l5l |. 

The paper is organized as follows. In the next section 
we repeat the basic equations for the evolution of ten- 
sor perturbations (gravitons) and we explain why it is 
not straight forward to include the velocity term of the 
boundary condition. In Section IIIII we present the new 
approach and obtain the modified perturbation equations 
via a coordinate transformation which is such that the 
velocity term in the boundary condition disappears. We 
then quantize the system in the new coordinates. In Sec- 
tion llVl wc show numerical results for graviton production 
at relatively low velocities. In Section [V] we conclude. 
Technical details are deferred to appendices. 



II. A MOVING BRANE IN AdS 5 
A. The background 

In Poincare coordinates (x A ) = (i,x, y) with x = 
(x 1 , x 2 , x 3 ) and A = 0,...,4, the AdS 5 (bulk) metric is 



given by 
ds 2 = g AB dx A dx l 



L 2 
V 2 



[-dt 2 



5ijdx l dx 3 +dy 2 ] , (1) 



where i,j = 1,2,3 and L is the AdSs curvature radius 
which is related to the bulk cosmological constant by the 
5d Einstein equation, —A = 6/L 2 . The physical brane 
representing our (spatially flat) Universe is located at 
some time dependent position y = yt(t) in the bulk, 
and the metric induced on the brane is the Friedman- 
Robertson- Walker metric 



ds 2 = a 2 (r]) [-drf + S.^dx'dx 3 ] , (2) 

with scale factor a{rf) which is given by the brane posi- 
tion, 

L 



a{rj) 



Vb(t) 



(3) 



The conformal time r\ of an observer on the brane, is 
related to the bulk time t via 



di] = VT - v 2 dt = j^ 1 dt . 



(4) 



(5) 



Here we have introduced the brane velocity 

dy b LH 1 

v — —— = = and 7 = - 

dt VI + L 2 H 2 VT^v 2 

H is the usual Hubble parameter, 

H = = a- l H = -L-^v . 
a 1 dr] 

Its dynamics, as a result of the second junction condition, 
is determined by the modified Friedmann equation [l| 



(0) 



H 



3 



—) 



(7) 



where T is the brane tension, p the energy density on the 
brane, and we assume the RS fine tuning condition Q 



12 



3 

T 2 



Furthermore (see Q), 



K4 



IttGa 



6 



(8) 



(9) 



L 



pi 



(10) 



We define the string and Planck scales by 
1 r3 1 

Note that the RS fine-tuning condition is equivalent to 

(11) 



T T 2 
k 5 — k 4 Ij or — — — — . 

L Li 



Identifying K5 with the string scale is based on the as- 
sumption that this phcnomcnological model comes from 
string theory with one large extra-dimension L, the y di- 
rection, while all the other extra-dimensions remain of 
the order of the string scale, L s . In this case the 4d 
observed Planck scale is related to the string scale by 

Eq. cm. 
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B. The setup 

We consider a radiation dominated brane which moves 
frome the Cauchy horizon, y ~ 0, at t = — ootoa position 
2/6 (0) < Vs at t = 0, where it bounces and changes its 
direction. In a radiation dominated universe p oc a -4 oc 
y b (t) i . Defining 



r(t) 



(12) 



we have i?L = y / r(l 
yields 



r/4). Inserting this in Eq. (|5|) 



y b (t) = v(t) = ± 



Vr(i)(l + r(t)/4) 



1 



-(i)/2 



(13) 



Here the upper sign is chosen for negative times, when y 
is growing and the universe is contracting while the lower 
sign corresponds to positive times (expanding universe). 
At the bounce the maximal velocity, v(0) is reached cor- 
responding to the maximal radiation density given by 



r(0) 




v 2 (0) 



1-U 2 (0) 



(14) 



At t ^ the radiation density is 



r(t) = r(0) 



Vb{t) 

y b (0) 



Note that since the differential equation Jl3|) is first order, 
only one initial condition, e.g v(0) can be chosen. yb(0) 
is then determined by the implicit equation 



y(o) 



v(t)dt 



(Implicit because it contains y(0) also in the integrand.) 
Because of this complication it is simpler to choose the 
initial conditions at some early time, t- m < so that 
K^in) ^ 1- For t<ti n we can then approximate Eq. CG 
to 



y b (t) = y/r(i) = VKfeO 



y 2 b (t in ) 



with solution 



Vb(t) 



t 



t < U r 



(15) 



The initial condition (t{ a , r(ii n )) determines the bounce 
velocity v(0). 

In this first paper, where we mainly want to present 
the method how to transform mixed boundary condi- 
tions into Neumann boundary conditions, we simplify the 
background evolution by assuming HL C 1 or, cquiva- 
lently, r(t) <C 1 at all times. This is of course not a good 
approximation if the bounce velocity is high and we shall 



treat the brane motion correctly in Ref. [jj 
(fT3"| reduces to 

Vb{t) = VW)^y 2 b (t) 

at all times. The expression 

L 2 



With this 



(16) 



Vb(t) 



\t\+t b 



(17) 



with parameter % solves Eq. p6|) for all t ^ 0. Fur- 
thermore, it has the correct asymptotics and the bounce 
velocity is given by 



V (°) = 72= —PT- = V b 
l b 



L 2 



The kink at t — can be regularized by replacing \t\ by 
y/t 2 + 1 2 , where t c is a small regularization parameter. 
For \t\ -C t c this does not affect the dynamics but for 
\t\ < t c the velocity is reduced and it actually passes 
through zero at t = 0. For graviton frequencies with 
Lut c <^ 1 the particle production obtained is independent 
of this regularization. 



C. Tensor perturbations 

Allowing for tensor perturbations hij(t,x,y) of the 
spatial three-dimensional geometry at fixed y, the per- 
turbed bulk metric reads 



d^ 



y 2 



-At 2 + (Sij + 2h ij )dx i dx j + dy 2 



(18) 



Tensor modes satisfy the traccless and transverse condi- 
tions, h\ = dihj =0. These conditions imply that hij has 
only two independent degrees of freedom, the two polar- 
ization states • = x,+. We decompose into spatial 
Fourier modes, 



hij(t,x,y) = 



d 3 k 
(2tt) 3 / 2 £ 



E 



> e?,(k)/i.(*,l«k) 



+,x 



(19) 

where (k) are unitary constant transverse-traceless po- 
larization tensors which form a basis of the two polariza- 
tion states • = x,+. Since the problem at hand obeys 
parity symmetry, we shall neglect in the following the 
distinction between the two graviton polarizations and 
consider only one of them. We then have to multiply the 
final results for e.g. particle number or energy density by 
a factor of two to account for both polarizations. 

Here we only consider 4d gravitational waves. The 
5d metric has in principle five different spin-2 polariza- 
tions. Two of them are the ones discussed here. In ad- 
dition there are the two helicities of the so-called gravi- 
vector and a gravi-scalar (see, e.g. [IB|). The gravi- vector 
and the gravi-scalar obey exactly the same propagation 
equation as the 4d gravitational waves in the bulk, only 
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their boundary conditions arc different. In principle they 
would add to the results obtained here. In this sense our 
results are conservative, but since the different polariza- 
tion states do not interact at the linear level they can 
be calculated independently. These polarizations are ex- 
pected to contribute on the same level as the two consid- 
ered here. 

The perturbed Einstein equations and the second junc- 
tion condition lead to the boundary value problem 



dl + -& 



h(t, y; k) = in the bulk (20) 



and 



7 (vd t + d y ) h\ 



Vb[t) 







(21) 



describing the time-evolution of the tensor perturbations 
as the brane moves through the bulk. We introduce also 
a second, static brane at position y s , which requires the 
additional boundary condition 



dyhL = 



(22) 



Eq. (|20|) is the Klein-Gordon equation for a minimally 
coupled massless mode in AdSs , i.e. the operator acting 
on h is just the Klein-Gordon operator 



□ 



f 



-9 



d A [a 



-99 



AB 



(23) 



Equation (|2ip is a time-dependent boundary condition 
(BC) coming from the fact that the moving brane acts 
like a "moving mirror" for the gravitational perturba- 
tions. Only in the rest-frame of the brane do we have 
pure Neumann BCs. In a generic frame we have the 
Lorentz transformed BC which contains a velocity term 
vd t . 

We assume that the brane is filled with a perfect fluid 
such that there are no anisotropic stress perturbations 
in the brane energy momentum tensor, i.e. there is no 
coupling of gravitational waves to matter. If this were the 
case, the r.h.s. of Eq. ([21]) would not be zero but a term 
coupling hij to the matter on the brane, see Eq. (2.25) 
of [12J, would be present. 

The analogy to a moving mirror is actually not just 
a pictorial one. Transverse-magnetic modes of the elec- 
tromagnetic field in an ideal, i.e. perfectly conducting, 
dynamical cavity are subject to the very same boundary 
condition, see, e.g., [l7j . In this context, the boundary 
condition (|21|) is sometimes referred to as "generalized 
Neumann" boundary condition, a terminology which we 
also adopt here. If the cavity is non-perfect, then also 
in the case of the electromagnetic field, the right hand 
side of the boundary condition contains a term describ- 
ing the interaction of the photon field with cavity mate- 
rial, similar to the anisotropic stress perturbations for the 
gravitational case considered here. This suggests that a 
brane with no anisotropic stresses could be termed " ideal 
brane" . 



For the tensor perturbations the gravitational action 
up to second order in the perturbations reads 

L 3 



S h = 4. 



2k 5 



dt 



y„(t) V L 



\d t h\ 2 -\dyh\ 2 



-em 



(24) 



One factor of two in the action is due to Z2 symmetry 
while a second factor comes from the two polarizations. 
As we have shown in [T^], the BCs (|2 1I22[) are indeed 
the only ones for which SSh = leads to the free wave 
equation (f2TJ|) . (In principle also Dirichlet BCs, i.e. h 
vanishing identically on the brane, lead to a wave equa- 
tion in the bulk. But besides leaving no room for a non- 
trivial dynamics of the gravitational waves on the brane, 
these are not obtained from the Einstein equations in the 
bulk.) 



D. Dynamical Casimir effect approach 

The wave equation ([2"0"[) itself is not time dependent 
and simply describes the propagation of free modes. It 
is the time dependence of the BC (|2"Tj) that sources the 
non-trivial time-evolution of the perturbations. As it is 
well known, such a system of a wave equation and a time- 
dependent BC leads, within a quantum mechanical for- 
mulation, to particle production from vacuum fluctua- 
tions. In the context of the photon field perturbed by 
a moving mirror this goes under the name "dynamical 
Casimir effect" @,@|. 

In (Til . [l2| we have extended a formalism which 
has been successfully employed for the numerical in- 
vestigation of photon production in dynamical cavities 
fl8l [la. f20l| to the RS braneworld scenario. We have stud- 
ied graviton production by a moving brane, which we call 
dynamical Casimir effect for gravitons, for a bouncing 
braneworld scenario. 

However, in order to solve the problem, we have ne- 
glected the velocity term in the BC (pH]) . The ansatz 



h = a a (t)t 



t <j> a (t,y)+h.c. , Lu 2 ak = k 2 + m 2 a {t) 



a=0 



then leads to a Sturm-Liouville problem for the instanta- 
neous eigenfunctions (j) a (t, y) consisting of the differential 
equation 



<j, a {t,y)=ml{t)<]> a {t,y) (25) 



and Neumann BCs at both branes. The solutions of ([25 
respecting Neumann BCs at both branes are 

y s yb(t) 



Mt) = 



(26) 



M*,V) = N n {t)y 2 C 2 (m n (t),y b {t),y) 
with 

C v (m,x,y) = Y 1 (mx) J u {my)- Ji{mx)Y v {my) . (27) 
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They form a complete orthonormal system with respect 
to the inner product 

0/0 = 2/ -4Mt,v)Mt,y) =*a/j (28) 

and the completeness relation implies 



(29) 



The factor two accounts for the Z2 symmetry of the bulk. 

In (l2^ we call <po and <p n the zero-mode and Kaluza- 
Klein (KK)- mode solution, respectively. Here 4>o is the 
massless mode, m = 0, which reduces to the usual 3 + 1 
- dimensional graviton on the branc. The KK masses 
m n 7^ are determined by the BC at the static branc, 
see, e.g. [H, HlJ for more details. 

Due to the completeness and ortho-normality of the 
functions {4>a} at any instant in time, any general solu- 
tion of l[2"5]) subject to Neumann BC's can be expanded 
in these instantaneous eigenfunctions. If we add the term 
vdt to the boundary condition this feature is lost and we 
can no longer expect to find a complete set of instanta- 
neous eigenfunctions. 

However, since the entire effect disappears when the 
velocity tends to zero, neglecting a term which is first or- 
der in the velocity seems not to be a consistent approach. 
This problem prompted us to search for another descrip- 
tion allowing us to treat the boundary condition (|2ip in 
full. 



III. GRAVITON PRODUCTION IN A 
TIME-DEPENDENT BULK WITH A MOVING 
BRANE 

In this section we introduce a new time coordinate 
which is chosen such that the velocity term in the bound- 
ary condition disappears but the mode equation for 
the instantaneous eigenfunctions still remains the Bcsscl 
equation (|2"B")) with its solution given by Eqs. (|2"6")l and 
([27f . We then extend the formalism of [l2| to this case 
and shall see that for small velocities our previous results 
are not modified. 



A new time coordinate 



We introduce new variables (x A ) = (r, x, z) given by 



r(t,y) = t + s(t,y) , z = y . 



(30) 



The idea is to find a function s(i, y) such that r — ► t for 
all y, when v — ► and that the junction condition (|2ip 
reduces to a normal Neumann BC in the new variables. 
We can then use the mode functions (|26|) and ([27)) to for- 
mulate the problem quantum mechanically One might 
first be tempted to make a y-dependent Lorentz trans- 
formation to the rest frame of the moving brane, but 



actually this does not lead to Neumann BC's in our case 
as the transformation induces new terms in the metric. 
We therefore first leave the function s(t,y) completely 
general and formulate the conditions which have to be 
satisfied in order for the new BC's to be purely Neu- 
mann. 

In (t, x, z)-coordinates, the brane trajectory is given 
by the implicit equation 



Zb{r) = y b [t(r, z&(t))] 



(31) 



Once we have specified the function s(t, y), the new brane 
trajectory Zb{r) can be found. This is done numerically 
since neither s(t, y) nor the inverse i(r, z) of (|30|) exist in 
closed form. As in fl2l ] we restrict ourselves to brane mo- 
tions where asymptotically, i.e. for t — * ±00, the physical 
brane approaches the Cauchy horizon (j/j — ► 0), moving 
very slowly (v — ► 0). 

The new metric given by 



Z.A-B 



ds =gAg(T,z)d£ x 



(32) 



is time dependent and contains non- vanishing cross terms 
go z . The explicit expression is given in (|A3[) . We now 
show that the function s(t,y) can be chosen such that 
the time-derivative term in the boundary condition dis- 
appears. 

In the coordinates defined in Eq. (|3H|) . the junction 
condition (f2"Tj) becomes 



v(t) 



dr dr 
dt dy 
v(t){l + d t s(t,y)} 
+d z h(r 7 z) 



d T h(r, z) + d z h(r, z) = 

d y s(t, y)] d T h(r, z) 
- at z = Zb(r) . 



(33) 



In order to obtain Neumann boundary conditions, we 
require that the term in square brackets vanish at z&(t). 
This leads to the condition 



d y s(t,y) 



l + d t s(t,y) 



v(t) 



(34) 



V=Vb(t) 



for the function s(t,y). Furthermore, we want to main- 
tain the Neumann BC at the static brane y s = z s . This 
yields the additional condition for the function s(t,y) 



d y s(t,y)\ y= y s = 



(35) 



Hence, if we can find a function s(t, y) which satisfies 
Eqs. ([34]) and ([35]) . the junction conditions in the new 
coordinates reduce to Neumann BC's 



d z h(r, z) = at z = z s and z = 2&(r) 



(36) 



To find a suitable function s(t,y) we choose the sepa- 
ration ansatz 



leading to 



v(t) 



s(t,y) = f(t)a(y) 

f(t)d y a(y b (t)) 
l+(d t f(t))a(y b (t)) 



(37) 



(38) 



i 



For the transformation ([30)) to be regular, we have to 
require 1 + d t s(t,y) = 1 + f t {t)a{y) £ V(t,y). If we 
choose a such that d y a(yi,(t)) is bounded from below, 
< A < d y a(yb(t)), this ansatz ensures the required 
asymptotics, f(t) — ► for v(t) — ► 0. In addition we need 



d v cr{y)\y=y B = 



(39) 



The function f(t) is determined by the differential equa- 
tion 



df(t) 1 <r'(y b (t)) 
dt v(t) a(y b {t)) 

A simple choice for a(y) is 



a{y) = l + — I 1 



/(*) 



1 



*(!&(*)) 



= 



(40) 



, do = const., (To > 1 , (41) 



so that 1 < cr(y) < 2. With this, condition (|35|) is auto- 
matically satisfied. In addition, we want the branc colli- 
sion, i.e. the bounce to happen at the fixed time r = 0. 
For this we chose the initial condition 



(42) 



/(t = 0) = 0. 



Since f(t) — > for v(t) — > 0, the transformation (|30l 
satisfies 



t — > i for i — > ±oo and r(t = 0,y) = 0. 



(43) 



For the first of these equations we use that v(t) — > for 
i — > ±oo and the form of the differential equation (|40|) . 

With (f3"5| . /(0) = implies that the velocity vanishes 
at the brane, u(0_) = «b = u(0+) = 0. Hence the ve- 
locity does not jump from a large value v b to — v b at the 
bounce but it evolves very rapidly but smoothly from 
a high positive value v max = v(—e) to a large negative 
value —Umax = u(e), e > and small (see Fig. [5]), like the 
regularized branc motion proposed in Section III Bl Con- 
firming that the results are independent of the choice of 
do is of course a crucial test. 

The coordinate transformation maps the problem of a 
moving brane in a static bulk {l} onto the problem of a 
brane moving according to (|3ip in a time-dependent bulk. 
At first glance a further complication of the problem. Its 
benefits, however, will become clear in the next sections. 
The transformation of the metric is given cxplicitely in 
Appendix A. 



B. Wave equation 

Transforming the Klein-Gordon operator (|2"3")l to the 
new coordinates x A , we obtain the wave equation 



,9i (r, z)d 2 T + g 2 (t, z)d T - 2s 2 (r, z)d z d T 



+ -d z - dl 



h(r, z) = 



(44) 



FIG. 2: The velocity in the new coordinates, 



Note that 



this is a coordinate velocity, not a physical quantity. It is easy 
to check that qab dXb d ^ < at all times, hence the 

physical velocity remains timelike also if becomes larger 
than 1. The maximal velocity for this case is Umax = 0.3. 



The definitions of the functions gi(r,z), giir^z) and 
S2(t, z) in terms of the coordinate transformation s(t, y) 
are given in Appendix [A] These functions manifest that 
the bulk itself is now time-dependent and that the met- 
ric is no longer diagonal. In Poincare coordinates the 
non-triviality of the time-evolution of the perturbations is 
purely a consequence of the time-dependent junction con- 
dition, no time-dependent functions enter the wave equa- 
tion (I20p . Our coordinate transformation which trans- 
forms the generalized Neumann BC into a pure Neu- 
mann BC, induces explicite time-dependence in the wave 
equation itself. What is important, however, is that in 
(|44[) . in the instantaneous rest frame where we neglect 
time derivatives, we get the operator (|25[) of the original 
Bessel equation with normalized solutions an d (|27[) . 
We just have to replace the variables (t,y) by (r, z). 

Writing the action (|24|) in terms of the new coordinates 
yields 



S = 4 / dr 



El 

2k 5 



d A k 



Zs dz 1 



26 (T) 



Z A 1 + Si 



(45) 



gi\d T h\ 2 - 2s 2 Re [(d T h)(d z h*)] - \d z h\ 2 + k 2 \W 2 



Using the expressions for s\, S2, gi and g2 given in Ap- 
pendix [A] it is readily shown that the variation of (|45[) . 
demanding Neumann boundary conditions at the brane 
positions, leads to the wave equation (|44|) . 

In the next subsections we take the action (|4"5"|) as the 
starting point to set up the dynamical Casimir effect for- 
mulation of graviton production along the same lines as 
in [ll|, [l2j . For the physical interpretation of gravitons 
we are using the fact that asymptotically, when the ve- 
locity of the brane goes to zero, the action ([45]) and the 
wave equation (|4"4"]) reduce to (|2"4")) and (|2T)|) , respectively. 
However, in the new coordinates, the junction conditions 
are always simple Neumann boundary conditions. 
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C. Mode decomposition and Hamiltonian 

As a basis for a mode decomposition we chose the 
cigcnfunctions {0 a (r, z)} obtained by replacing (t,y) — > 
(t,z) in (gni) and P?]), As in [T3] we call O and 4 the 
zero-mode and KK mode solution, respectively. For a 
brane at rest, and hence r = t, the solutions <pQ and 
(^>j do indeed represent the physical zero mode and the 
KK modes, see, e.g. 0]. When the brane is moving, 
however, these solutions are 'instantaneous modes', pro- 
vided that the boundary condition is Neumann. This 
approach is widely employed in the context of the dynam- 
ical Casimir effect, see [IH, [ljlHtl and references therein. 
Here, working in the (r, z)-coordinates, the modes (|26|) 
and (|2"T|) . are proper eigenfunctions respecting the full 
junction condition which we have reduced to a Neumann 
BC. At early and late times, i.e. asymptotically \t\ — > oo, 
where the brane velocity tends to zero, these eigenfunc- 
tions agree with the physical cigcnfunctions correspond- 
ing to the zero mode and the KK modes. Since the eigen- 
functions {0 Q (r, z)} form a complete and orthonormal 
set. and satisfy the correct junction conditions at both 
branes, we may decompose the graviton field in ^> a 's and 
the pre- factors q a ,k{ T ) become canonical variables which 
can then be quantized (IT 



h(r,z,k) 




3a,k(r)</> a (r, z) 



(46) 



Our coordinate transformation and the expansion (|46p 
satisfy two major requirements. First, the expansion 
(|4"6")l is consistent with the full junction condition (gen- 
eralized Neumann BC). This overcomes the problem of 
our approach in [ll], [lj] • Secondly, even if at arbitrary 
times the q a ,kS cannot a priori be identified with phys- 
ical modes, asymptotically, i.e. when the brane moves 
very slowly, they do represent the independent physical 
graviton modes. This allows us to introduce a proper no- 
tion of particles and vacuum states for asymptotic times. 
Initial and final vacuum states are then linked by the 
time-evolution of the q a ,k§ exactly as in fl2T |. 

We divide the wave equation ((44)) by <?i in order to iso- 
late the second time derivative and insert the expansion 
(|4"6"|) . Note that g\ — > 1 for \t\ — > oo and for a sufficiently 
large choice of ctq, g\ > at all times. As we shall see 
below, this is also needed for the Hamiltonian to be pos- 
itive at all times. Inserting the expansion (|46[) into (|4"4")) . 
multiplying it by <pp and integrating over 2 f**,^dz/z 3 
leads to a system of differential equations for the q a ,k 
which has the same form as the one of Refs. (Til. [l2j. 

Qa,k(r) + ^ [^/3a( T )9/3,k(T) + Bp a {r)qp^.{T)) = . 



(47) 

The explicite expressions for the time-dependent cou- 
pling matrices Ap a and Bp a are given by integrals over 
the bulk which are rather cumbersome. The details 
can be found in Appendix [B] Inserting the expansion 



(|46| into the action (|45|) we obtain the Lagrangian L(t) 
in terms of the variables q a ,k(T~). We can then define 
the canonical momenta p a ,k = dL/dq a ^ from which, 
by means of a Legendre transformation, we derive the 
Hamiltonian 



J a/3 



pP/3,-k 



+q a ,k 



-(M 0a - Sp a ) [qf3,kPa,k + p 



qp, 



(48) 



The matrices E~p, V a p, Mp a and Sp a are given ex- 
plicitely in Appendix |Bj It is important to note that 
E~p is positive definite as long as g\ > and 1 + s\ > 0, 
which we have to require for our approach to be consis- 
tent. In the old treatment, E a p was the identity ma- 
trix and the couplings Sp a and V a p were missing. They 
are due to the time-dependence of the bulk spacetime 
in the new coordinates and therefore originate from the 
term vdt of the boundary condition in Poincare coordi- 
nates. The coupling matrix Mp a which is also present in 
our previous treatment comes from the time dependent 
Neumann BC. Finally, the time dependence of the bulk 
volume z s — Zj,(t), induces the time dependence in the 
frequency 0J a ,k (squeezing effect, see [H[). 

All the coupling matrices tend to zero when v — > 0. 
But we have not been able to show that the new 
couplings, E~l — S a p, V a p and Sp a are parametrically 
smaller than Mp a , e.g. that they are of order v 2 . There- 
fore, the result that the par ticle production obtained in 
our previous treatment [ill . Il2| is not modified if the 
velocity is sufficiently low is not evident and has to be 
checked numerically. 



D. Quantum Generation of Gravitons 

The quantization procedure goes along the same lines 
as in [121 ]. The canonical variables q a ,k{T~), Pa,k(j) 
and the Hamiltonian H(t) are promoted to operators 
qa,k{r)i Pa,k(T) and H(t), subject to the usual commu- 
tation relations. In the Heisenberg picture where the 
time evolution of an operator O is determined by 



0(T)=i[H(r),0(r)} + 



dO(r) 
dr 



cxpl. 



the operators q a ,k{T) andp Qj k(r) satisfy the same Hamil- 
tonian equations of motion as their classical counterparts, 
i.e. 

dH dH 



qa.k 



dp a ,: 



Pa,k 



4a, k 



(49) 



Remember that we assume that asymptotically, \t\ — > oo, 
the brane is at rest, i.e. the brane velocity vanishes and 
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both coordinate systems agree, r = t. Wc extend this 
notion of asymptotic behavior by introducing two times, 
ti n and tout , and we shall assume that the brane is at rest 
for t < tin and t > t ou t, respectively. This corresponds 
to a scenario where the motion of the brane is switched 
on and off at finite times. Such a brane dynamics may 
seem rather artificial from a physical point of view, but 
what is important for us is that before ti n and after t out 
no significant particle creation takes place. Numerically, 
we test this by varying t m and i out and choosing them 
large enough so that the particle number is independent 
of the value chosen. 

In the (r, z)-coordinates, the brane is then at rest for 
times t = t < = tm and r = t > r out = t out , 
respectively. When the brane velocity is zero, the ma- 
trix E a fj{T) defined in Appendix [Bl becomes the identity, 
E a p{T) ^| r |^oo ^a/3j and all other matrices which rep- 
resent the coupling terms vanish identically in this limit. 
Consequently, for asymptotic times the Hamiltonian re- 
duces to the familiar form of a collection of independent 
harmonic oscillators, 



H' 



with 



1 / out _ 

~~ 2 



(50) 



(51) 



Pa,k = 9a, -k • 

We have introduced the notation 

0J™ k EE LO a A T < Tin) , <^ = CJ a ,k{ T > T °ut) ■ (52) 

Following [l2| , we decompose q a ^ in creation and anni- 
hilation operators, 



g Q ,k(r) = ^2 

p J2uj 



a l3,k e a ,k\ T ) + a /3,-k e a,fc V T > 



inf (0)*, 



:> k 



(53) 

which are defined via a™ k |0,in) =0 V a,k. The 
initial vacuum state |0, in) is the ground state of the 
Hamiltonian (|50[) for times r < Tj n . This initial state is 
linked to the final vacuum state defined by a° u £|0, out) = 
; V a, k, by means of a Bogoliubov transformation (see 

m) 

= E [A,/3 I fc(r »t)S? lk + B^AroutKU] (54) 



IV. NUMERICAL RESULTS 

In order to solve the equations of motion l|47j) numeri- 
cally, we transform the system to a first order system and 
introduce a mass cutoff, n max , i.e. we neglect all modes 
with masses higher than m„ max , in other words q a = for 
ot > fimaxi along the lines explained in detail in Ref. [13]. 
Modes close to this cutoff are of course seriously affected 
by it as is seen in Figs. 2] and [5] We have tested the 
stability of the results for modes n <C n max by varying 
the cutoff. The stability of the zero mode is illustrated 
in the lower panel of Figs. [3j and [4j Typically modes 
with n <,0.7n max can be trusted. An indication for this 
is also the Bogoliubov test shown in Fig. [6] and discussed 
in Ref. [13]. 

The first order system is given explicitly in Ap- 
pendix IB 31 and differs from the original one in [l2j only 
by additional mode couplings. We have compared our 
new results with those of Ref. [TJ [13] and find excellent 
agreement at low velocity, t> m ax~0.1. This is illustrated 
in Figs. [3] and [4] At bounce velocities v m ax^0.5 we do 
find differences as expected, but these results cannot be 
taken literally since for these velocities the low energy 
evolution of the scale factor adopted in this work is no 
longer sufficient. We will present the full high velocity 
results in a forthcoming paper fTHl ]. 

The agreement between the old results for pure Neu- 
mann boundary conditions and the new ones with the 
generalized Neumann boundary conditions is similar for 
other values of y s and k. 

In Fig. [5] we show the KK spectra for v ma x = 0.1 and 
» max = 0.3 for the wave number k = 0.1 and position of 
the static brane, y s = 10. In this case, the analytic ap- 
proximation derived in Ref. [13] which is valid for m n < 1 
can only be trusted for the first two modes. The slight 
difference between the old and the new spectra towards 
the end, i.e. for m n > 10, is due to changes how wc 
numerically evolve the solutions through the bounce (see 
Appendix C). This affects the sensitivity of the solutions 
to the cut-off. What we observe here as a slight bending 
of the spectrum for generalized Neumann BC's is also 
found in our previous approach if we increase the num- 
ber of modes: compare to the n max = 100 results shown 
in Fig. 15 of [l2| and the discussion related to Fig. 25 of 
[l3]. The drop in the final part of th e sp ectrum is just 
an artefact of the finite cut-off (see [l^] for a detailed 
discussion) . 



which determines the number of produced gravitons (for 
each polarization) 



N°S=J2\ B ^MTout)f 



(55) 



As we have discussed in detail in [121 ] , the graviton num- 
ber after the time r out , ([55]) . represents a physically 
meaningful quantity . 



V. CONCLUSIONS 

In this paper we have derived a method to calculate 
graviton production in bouncing AdSs brancworlds by 
the dynamical Casimir effect taking into account the 
full generalized Neumann boundary condition. We have 
achieved this by transforming to a new time coordinate, 
in which the generalized Neumann BC become ordinary 
Neumann BC. Wc have shown numerically that for low 
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FIG. 3: The final graviton spectrum for three-momentum 
k — 0.01, brane separation y 3 = 100 and bounce velocity 
i'max = 0.1. The top panel shows the final KK mode spec- 
trum and the lower panel depicts the time evolution of the 
zero mode. What we plot here is a kind of instantaneous 
particle number (see Appendix C of [U)- The numerical 
result for the KK spectrum is compared with the old one 
(shown in Fig. 13 of [ijj )• Like there, lengths are in units of 
L and momenta/masses in the units of L~ . For low veloc- 
ities Umax < 0.1 the new spectra (generalized Neumann BC) 
are identical with the old ones (Neumann BC) within the nu- 
merical error which are estimated by the Bogoliubov test (see 
Appendix). No.k(i~) is shown for two cut-off parameters n ma x 
to underline stability of the solution. 



bounce velocities, w max 0.1, the number of generated 
particles agrees with the one from the simpler treatment 
which neglects the velocity term in the boundary con- 
dition. Since this term is of first order in the velocity, 
we believe that our result is not obvious. Furthermore, 
the method developed in this work can be used to cal- 
culate particle creation for branes moving at arbitrarily 
high velocities. In this case, one will have to take into 
account the modification of the Friedmann equation at 
high energy, HL >, 1 . This is the goal of a forthcoming 
paper [l5j . 

In this work we have not derived new physical results, 



FIG. 4: As Figure [3] but for i) max = 0.5. For this velocity 
we do see a difference of about 10% between the previous, 
inconsistent approach and the new generalized Neumann BC 
for both, the KK modes as well as the zero mode. Again, the 
4d graviton number is shown for two cut-off parameters n max 
to indicate numerical stability. . 



but we have developed a new, fully consistent method 
to calculate graviton production due to the motion of a 
braneworld. Our method overcomes a shortcoming of our 
previous treatment [ill Il2| , and we have verified that at 
low brane velocity, w max ;$0.3 the previous results are not 
affected. 
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FIG. 5: Final graviton spectra for three-momentum k = 0.1 
and brane separation y s = 10 for the bounce velocities 
v max =0.1 and 0.3. Again, the new (generalized Neumann 
BC) numerical results are compared with the ones of the pre- 
vious inconsistent approach (Neumann BC), see Fig. 15 of 
and the agreement is excellent in the regime m n < 13, 
where the numerics can be trusted. 



and its inverse is 



d(r,x,z) 



i+a t s u u u i+d t s 

10 

10 

1 

V o ooo i j 



Under this coordinate transformation the AdSs metric in 
Poincarc coordinates given in (fTJ) transforms to 



gAB<lx A dx 



((T-y 9 T-i) AB dx A dx E 
^ 1 < dr 2 



(A2) 



.(1 + Sl (r,z))2 . 
+2s2(t, z)drdz + gi(r, z)dz 2 ^j + 5ijdx l dx^ 

We introduce the functions 

Si (t,z) = {d t s)(t{T 7 z) 7 z) (A3) 

s 2 (r,z) = (d z s)\ t=const (t(T,z),z) (A4) 

s n (r,z) = (d 2 s)(t(T 7 z) 7 z) (A5) 

S22 (t,z) = (d 2 z s)\ t=const (t(T,z),z) (A6) 



and 



9i{t,z) = (l + si(r, z)) 2 - s 2 (t, z) 2 



(A7) 



92{t,z) = sn (r, z) - s 2 2 (t,z) + -s 2 (t,z) , (A8) 

of which gi and g 2 will be used in Appendix [Bl 
The determinant is 



'L\ W 1 

g = det{g AB ) = - ( - — — 

z J (1 + si(r,z)) 2 



(A9) 



APPENDIX A: THE COORDINATE 
TRANSFORMATION 



APPENDIX B: DETAILS ON EVOLUTION 
EQUATIONS 

1. Wave equation 



The Jacobian T of the transformation 



The coupling matrices which determine the mode evo- 
lution equation (|47[) are given in terms of the following 
bulk integrals: 



(t,y) h-> (t = t + s(t,y),z = y) 



reads 



8(t,: 



d(t,x,y) 



1 + d t s 











d y s 







1 

















1 



















1 



















1 


) 



A, 



a/3( 



dz 



2<j> a 



92(t,z) 



<f> a (r, z) 



2s 2 (t,z) 



dz 



9i{t, z) 

Mt,*) (Bl) 



z b{r) 



<Pa(T, Z) -\ r(pa{T, Z) 



<A1> 2 «<^')< (T ,, )+ -ii<4,„ (r , S ) 



9i(t,z) 



9i(t,z) 



9i(t,z) 

^(t,z) (B2) 
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with 



Ua,k( T ) = \Jm 2 a {T) + k 2 



(B3) 



The over-dot denotes the derivative w.r.t. the time r and 
a prime stands for the derivative w.r.t. the coordinate z. 
Compared to our former work [l2j , the present problem is 
more complicated due to the additional couplings which 
are caused by the time-dependence of the bulk space- 
time. Also the Lagrangian and Hamiltonian equations 
for q a j. are more complicated. Furthermore, the func- 
tions si,S2,gi and gi are only known numerically. This 
induces additional numerical difficulties. Note also that it 
is important that g\ does not pass through zero for these 
integrals to be well defined, hence gi(r, z) > Vr, z. This 
is, however, easily achieved with our ansatz (|37l |4T|) for 
s(t,y) if we choose ctq sufficiently large. 



2. Lagrangian and Hamiltonian formulation 

Inserting the expansion (|46|) into the action (|45f leads 
(for each of the polarizations) to the Lagrangian 



l{t) = \ J rf3fc 5Z [ E <*p<ictM0 

a/3 



-k 



+ {M a f3 - IC a p){q a ^qp-k + q a -kq0,k) 
+ (A/" Q) 3 — V a /3 — Qaf3 — w a/3,fc)9a,k9/3,-l 



containing several time-dependent coupling terms. In de- 
tail, these read 



E a p{r) 

/C a| g(r) 
Q a p{r) 



r 


dz 


9i{t,z) 


b a (r, z)(/)p(T, z) 




z 3 


1 + si(t, z)' 


r 


dz 


9i{t,z) 


b a {r, z)4>p(T, z) 




z 3 


i + si(r, z y 


1" 


dz 


9i{t,z) 


j>a(T, z)<fip(T, Z) 




z 3 


i + si(t, z y 


r 


dz 


«2(t, z) 


t>' a (T,z)(l)0(T,z) 




z 3 


1 + si(r, z) 


r 


dz 


S2(T,Z) 


_0 Q (r, z)</^(t,z) 


)z b (r) 


z 3 


l + «i(r, 2) 



+^(t,«)^(t,z)] 
dz s[(t, z) 



*»(r) ^ (! +Sl(T,z)) 

+0' q (t, z)^(t,z)] 

i(m2( T )+m2(r))+fc 2 
dz <t> a {r, z)4>p(t,z) 



2 [0 Q (r, z)^(r, z) 



zb(r) 2:3 l + si(r,z) 



Since we require <7i(t, z) > and 1 + Si(r, z) > 0, the 
matrix E a p is positive definite. This is important for 



the above Lagrangian to lead to consistent second order 
equations of motion for the variables q a ^ (no ghosts). 

The equation of motion for the canonical variables ob- 
tained from the Euler-Lagrange equations become 

[j5 a7 g a ,k + E ai q a ^ + (M - ICj q a ^ 



+ [{M-JC) ai -(M-K,) ia 

(J\faj 7^a7 Qaj ^07. A:) Qa,k 



<7a,k 

= . 



(B4) 



Note that all the matrices introduced above apart from 
E a p tend to zero when v — > 0, i.e. for |r| — > oo. In 
this limit £^ Qj g tends to 6 a/ 3 so that Eq. (|B4[) becomes 
the free, uncoupled mode evolution equation in this limit 
as is expected. Introducing the canonically conjugate 
variables 



Pa,k 



dL 
dq a ,k 



^ [E a pqp,--k + {Mf3a - fCi3a)qp--k] 



(B5) 

leads by means of a Legendre transformation to the 
Hamiltonian (j48|) with coupling matrices 



V a0 (r) 



dz 1 

z b (r) ^l + Si 



9i 



- SI (t> a (r,z)(l)g(T,z) 



fc 2 si(r, z)<j> a {r, z)4>p(t, z) 
dz 



z b {r) 



j0a(T, z)(j)f}(T, z) 



S a f3(T) 

Thereby E~l is the inverse of E a p, i.c 



z„(r) Z A 31 



(B6) 
(B7) 

(B8) 



E 



a(3 



' s dz 1 + si(t, z) 
(r) ^ 3 9i(r,z) 



(f> a (T,z)M^z) . (B9) 



The Hamilton equations 
dH 

Qa, k 



, Pa,k 



an 

<9<7a,k 



(BIO) 



then provide the equations of motion for the variables 
Q a ,k and p aM . 

Using certain relations of the coupling matrices follow- 
ing from the completeness (|29| and ortho-normality ([28]) 
of the functions <fi a and the properties of the functions 
si, S2, sn and S22 one can show that the three systems of 
equations (|4"T|) . (|B4[) and the Hamilton equations (|B10|) 
are consistent with each other, i.e. one system follows 
from the other one. This seems to be at first sight a 
rather trivial statement but we have to remind the reader 
that this is not the case in our previous work [ll], [T3| as 
we have discussed in detail in Section II. D of Ref. [12 1 . 
The new coupling matrices V a p, S a p and E a p — 5 a p are 
missing in our previous work. Even though they do be- 
come very small when the brane velocity becomes small, 
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it is not evident that these new terms must be smaller 
than e.g. M a p, which also tends to zero with v. In other 
words, it is not straight forward to show that these con- 
tributions are, e.g., of order v 2 . 



3. Bogoliubov coefficients 

Performing the quantization as in Ref. [l2T ] we again 
transform to a first order system of equation. In the new 
coordinates the system of equations (3.34), (3.35) of [HJ 
is replaced by 

&(r) = e{-K^)+ c ^H4 7) w 



+ [^p(r)-cUr)]v { ;\r)} (Bll) 
■[^(r)+^(r)]^(r)} (B12) 



where 



i /i 



(B13) 



1 

2 

±- 



Mp a {r) - Spa (r) 



(M Q/3 (r) - S aP (r)) 



(B14) 



(Note that in Ref. [12| a factor of two is missing in the 
expression for Mf, in Eq. (B8), a simple misprint.) 



APPENDIX C: NUMERICS 

To compute the graviton spectra we have adapted the 
code described in Ref. [ljj to the new problem. Apart 
from calculating the new coupling matrices we also have 
to solve the differential equation (|40|) numerically to cal- 
culate the coordinate transformation and its inverse in 
order to determine Zb(r) via the implicit equation (|31|) 
as well as the functions si(r, z), S2(t, z), gi(r,z) and 
92 (t, z) which enter the integrals for the coupling ma- 
trices. For numerical purposes we have smoothed the 
function yb(t). Due to this implicit nature of the coor- 
dinate transformation, the calculation of these coupling 
matrices is numerically significantly more involved than 
in our previous approach. 



As in [12|, splines are used to interpolate the various 
matrix elements between time steps. The time steps used 
to produce the splines are not uniformly distributed but 
carefully selected to take into account the steepness of the 
time dependence of the couplings. Close to the bounce we 
use very short time steps to produce the splines (~ 10~ 6 ) 
while far away we can increase the step up to 0.2 (in units 
of L). Furthermore, due to the complex time dependence 
of some of the couplings very close to the bounce, ex- 
act integration of the matrix elements when propagating 
the solutions through the bounce is necessary in order 
to obtain satisfactory accuracy for large KK masses as 
in Fig. [5l In this way the bounce is taken into account 
as accurate as possible. This affects the speed of conver- 
gence of the solutions w.r.t. n max , leading to the behavior 
described below Fig. [5] 

As an indicator for the accuracy of our calculations 
we use the Bogoliubov test as described in Appendix D 
of (Eq. (D6)). This is presented in Fig. © for the 
I'max = 0.3 result given in Fig. [5l The quantity denoted 
by 'Bogoliubov test' and shown as solid line in Fig. [S] 
should ideally vanish. Given the complex nature of the 
numerical problem, the accuracy of the results is satis- 
factory for m n <, 10/L. 
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FIG. 6: Comparison of the final KK spectrum iV^.i and the 
corresponding quantity d n ,o.i(r ou t) given in Eq. (D6) of [l2l [. 
The quantity d„,o.i (r ou t) is supposed to vanish identically, 
see [12] • The comparison is a measure for the accuracy of the 
"max = 0.3 result depicted in Fig. [5] In the region of the spec- 
trum which is free from numerical artefacts, i.e. no depen- 
dence on cut-off n max , d n ,o.l (lout) is a t least about two orders 
of magnitude smaller than the physically relevant quantity 
indicating a satisfactory accuracy. 
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